x.c.b = cbind(exprSet.rorc,exprSet.rorb)
colnames(x.c.b) = c("RORC","Type","Gene","RORB","Type","Gene")
x.c.b = x.c.b[,c(1,4)]
library(ggplot2)
library(ggpubr)
p1 <- ggplot(data = x.c.b,mapping = aes(x = RORC, y = RORB)) +
  geom_point(colour = "red",size = 2) +
  geom_smooth(method = lm,colour = "blue",fill = "gray")
p1
p1 + stat_cor(method = "pearson",label.x = 0.5,label.y = 0.3)
